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Abstract This work is focused to study the development 
and control of the laminar vortex breakdown of a flow en- 
closed in a cylinder. We show that vortex breakdown can be 
controlled by the introduction of a small fixed rod in the axis 
of the cylinder. Our method to control the onset of vortex 
breakdown is simpler than those previously proposed, since 
it does not require any auxiliary device system. The effect 
of the fixed rods may be understood using a simple model 
based on the failure of the quasi-cylindrical approximation. 
We report experimental results of the critical Reynolds num- 
ber for the appearance of vortex breakdown for different ra- 
dius of the fixed rods and different aspect ratios of the sys- 
tem. Good agreement is found between the theoretical and 
experimental results. 
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1 Introduction 



The development of structural changes in vortical flows and, 



tensively investigated during tl 
li; 3 13; 14; 3 3 3D) 




particularly, vortex breakdown 
tensively investigated during the last years 

The characteristic and fun- 
damental signature of vortex breakdown is the appearance 
of a stagnation point followed by regions of reversed axial 
flows with a bubble structure when the swirl is sufficiently 
large. This structural change is also accompanied by a sud- 
den change of the size of the core and the appearance of 
disturbances downstream the enlargement of the core. Vor- 
tex breakdown is very important in several applications of 
Fluid Mechanics such as aerodynamics, combustion, nuclear 
fusion reactors or bioreactors. The presence of vortex break- 
down may be beneficial or detrimental, depending on each 

HQ 

particular application dllhll3h . 
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In order to explain the origin of vortex breakdown many Recently, different methods for controlling VB in closed 

proposals have been given. Some of them are based on in- flows have been proposed using different techniques such as 

stability mechanisms, while in other cases the proposed the- co-rotation and contour-rotation of the end-walls J13L the 

ories consider that hydrodynamical instabilities do not play addition of near-axis or at the end wall swirl llli . or tem- 

a significant role. In any case, it is widely accepted that in- perature gradients 1 lj]). Experiments show that such meth- 

stabilities arise after the occurrence of vortex breakdown. ods may increase or decrease the critical Reynolds number 



Vortex breakdown (VB) has been observed not only in 
open flows, but also in experiments performed in confined 
flows, for example, in closed cylinders (3). It is worth not- 
ing that the characteristics of VB in both cases are strongly 
similar, suggesting the possibility that the basic mechanism 
of VB is the same in both situations 

While considerable theoretical effort has been done in 



to develop VB. In all the mentioned cases, it is necessary 
to use auxiliary devices, such motors or heat sources. The 
aim of this work is to introduce a method simpler, from a 
practical point of view, than those previously mentioned for 
controlling vortex breakdown in closed flows. 

This paper is organized as follows. In Sec.|2]we analyze 
a simple model to predict the onset of vortex breakdown. In 
Sec. |3] we present the experimental observations. These re- 



, r-ijr. ■ n 1 i suits show that control of VB may be performed with the 

the study of VB emergence in open flows, less results have 

, 1 . - 1 n 1 1 n T .,. , presence of rods at rest, located along the axis of the cylin- 

been obtained tor closed flows. In a way, this may be ex- r e J 

. . , , - . . . . „ , . „ der, without the necessity of transferring swirl. The compar- 

plained by the fact that in the later case, the basic flow to J * 

, ,. , . , , 11 u ison between the experiments and the theoretical model is 

be studied is much more complex than those typically ob- v 



served in open ducts. On one hand, experimental measure- 
ments show that flows in open channels can be accurately 
described with the relatively simple q-vortex model o). On 
the other hand, an analytical description for closed flows 
which are considerable more complex than closed flows is 



given in Sec.|4] Finally, in Sec.|5]we present a summary and 
the conclusions. 



2 The model 



not known ( 3 



8|). However, from an experimental point of As already mentioned in the introduction, the similarity be- 



view, experiments performed with flows confined in cylin- tween vortex breakdown in open and closed flows is striking, 

ders are very attractive because they are simpler to control. In particular, experiments and simulations show that, the ve- 

Due to the great number of practical implications of VB, the locity fields around the point where VB develops, have the 

development of mechanisms for controlling its emergence is same structure in open and closed flows. The similarity of 

of considerable interest. both situations suggests that the formation of VB depends 
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mainly on the local velocity field, and the details of the flow This vortex has a core of radius cq which rotates like a solid 
far from the location of VB can be neglected. For example, body with constant angular velocity a and an outer irrota- 

in figure^ we can observe (a) the experimental streamlines tional region. The whole vortex has an axial motion with 

fl 

just before the appearance of VB from Ref. (3) and (b) a velocity Uq. In a cylindrical system of coordinates (r, d,z) 
sketch of an open flow with divergent section. In subsection with the z-axis in the direction of the duct the velocity field 
12.11 we present a simple model to predict the onset of VB is written as: 
in open flows. We shall employ this model to describe lo- 
cally the enclosed flow and to study the influence of local 
geometrical changes on VB. In particular, the effect on VB 
development of a fixed rod located at the axis of the cylinder 
is discussed in subsection !2.2l The results will be compared 
with the experimental observations in subsequent sections. where r . g ^ circulation in the outer region and in order to 



V r = 0, (1) 

or < r < cq, 



Vfl = < 



r/r c <r<R , 
V z = Uq. 



assure the continuity of the velocity at r = cq it must satisfy 

r = oc 2 . 



2. 1 Failure of quasi-cylindrical approximation and VB 

In the past, it has been shown that simple theoretical vortex 
models of open flows may predict behaviors very similar to 
those occurring in vortex breakdown. Various authors asso- 
ciated the emergence of VB with the failure of the quasi- 

II 

cylindrical approximation in inviscid theories This ap- 
proach gives a prediction about the position of VB that is in 

good agreement with experiments (4). A classical approach v z r ^ r ' 

1 dW 

based on an axisymmetric inviscid analysis was given by v r = — . (3) 

r az 



Now, in order to solve the fluid motion downstream, re- 
gion II in Fig.n we introduce the axisymmetric streamfunc- 
tion W related to the velocity field by means of 



1 3. 



Batchelor (20). In this model, that can be applied outside 
the boundary layer, it is assumed that the fluid in steady mo- 
tion passes through a transition from a duct of radius Rq to 



From Euler's equation the following equation for the stream- 

i3 



function is obtained I 



another duct of different radius R as it is shown in Fig. \T\b). d.ldW. d 2 W 2 dH T .dK 

r 2l(z^7) + ^T = r ^3777 ( 4 ) 



Based on Batchelor's model, we assume that upstream, 



dr K r dr ' dz 2 dW dW 



region I in Fig. Gib), the fluid moves inside a duct of radius where H = \v 2 + £ , K = rvg, p is the density, p is the 

z p 



Ro following a Rankine vortex with uniform axial velocity, pressure and v 2 is the magnitude of the velocity. For the flow 
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given by Eqs. ([0, we have that automatically satisfied because the flow has no radial com- 

ponent. Therefore, the downstream cylindrical flow in the 

2a 2 1 9 2a 

H = —W + -U Q 2 ,K=—W,forO<r<c (5) ^ u • t u A - f t u t 

Uq 2 Uo rotational region (r < c, where c is the radius or the rota- 

tional core) is given by 

Thus, Eq. (|4j takes the form 



1„. „. 4a 2 „, 2a 2 



%r - -% + W zz = - -^f + ^-r 2 (6) 

where the subscripts, r and z, denote derivatives. The rela- 
tions Q hold also downstream for the steady flow, so that 
Eq. determines the cylindrical flow downstream. The gen- 
eral solution for the downstream rotational region of the flow 
for any kind of cylindrical ducts is 

W( r )= 1 -Uor 2 +AF 1 (r r )+BY l (y r ) (7) /?j-c 2 „ 

2 Vz = R2- C 2 U °> (12) 

where F\ and Y\ are the Bessel functions of the first and sec- where, using mass conservation, 
ond kind respectively and y=2a/Uo (20). The downstream 

(cl-c 2 )U 

flow in the irrotational region is given by (l c ) 

and S is the swirl parameter, defined as S = Ry= 2 ^-. Us- 
v _ q ^ ing the continuity of the pressure, we obtain the following 

implicit equation that gives the value of the core radius c, 



V r 


= 0, 




ve 


= or- 


AS , ,s 


v z 


= t/cH 


AS S 
-^ M R 



(11) 

whereas in the irrotational region ( c < r < R), the flow is 
written as 

Vr = 0, 

ve = r/r, 



ve = r/r, (9) 
v z = U, (10) (R 2 Q -c 2 ) |(c 2 -c 2 )/ (|c 



(R 2 -c 2 ) IcJ^c) 



(14) 



where U is a constant that can be determined using conser- 
vation laws. We solved this equation for different values of the pa- 
Now, using the fact that downstream the fluid is con- rameters. When the value of S is increased, it may happen 
strained to move inside a cylindrical duct of radius R, in that two branches of solutions for c collide and disappear, 
order to avoid the divergence of W at r = it follows that In Fig.|2]the curve labelled (a) represents the ratio c/R as a 
B = 0. The no mass flow condition at the solid boundary is function of the swirl S. We observe in this figure that for each 
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curve (corresponding to the different values of d, parameter by the Eqs. Jilt . Thus the radius of the core c is to be de- 
to be introduced in the following subsection) there exists a termined with the Eq. i ll 4b . As already mentioned, for some 
critical value, above which the solution does not exist. range of values of Rq, Co and Uo, cylindrical solutions do not 

exist in region II if S is above a critical value S c (see Fig. 



13. According to Batchelor's criterion, this means that VB 
takes place inside region II, when the swirl parameter S is 
larger than a critical value S c . This behavior is in qualitative 
agreement with VB inside the closed cylinder. 

2.2 Effect of rods presence at the axis 



In Batchelor approach J2CJ) . the disappearance of these 
solutions is interpreted as the signal of VB emergence. We 
discuss now a way to support this interpretation. As already 
mentioned, the typical signature of VB is the formation of a 
stagnation point, followed by regions with reversion of the 
axial velocity. In the neighborhood of the stagnation point, 
a quasi-cylindrical description of the flow is doomed to fail- 
ure, because in this place the flow is strongly dependent on 

the axial coordinate z \§. This is revealed by the disappear- We now estimate how the changes in the geometry of the 
ance of cylindrical solutions. In the absence of VB a quasi- re g ion 11 affect the development of VB. More specifically, 
cylindrical description is suitable and, probably, there exists we consider the effect produced by the presence of cylinders 
a cylindrical solution to describe the flow. Then, it is rea- a l° n g the duct axis. Taking into account the local similarity 
sonable to assume that the critical value S VB for the emer- between open and closed flows, we extend the model intro- 
gence of VB and the critical value S c for the disappearance duced in the previous sub-section to this case. Using the for- 
of cylindrical solutions are very close to one another, i.e. malism described above, we obtain the following equation 
S VB ~ S c , and the identification of both may be a practical that determines the radius of the rotational core c in pres- 
criteria for VB. From now on, we use this criterion to es- ence of the inner cylinder of radius d 



{Rl-cD A U S S B n S S 

J o(^c)-—-Y {-c) = l, (15) 



with 



timate the critical conditions for the emergence of VB. We 

shall employ the above model of open flows to study qualita- (R 2 ~c 2 ) Uo R °^R 1 UqR°^R 
tively the flow in the closed cylinder. We shall assume that, if 
the open and closed flows are similar around the localization 

of VB, then the phenomena that takes place in both cases are , Uo r d i c l - c2 ) Y \{j d ) +cd 2 Y 1 (^c) 
similar. The flow inside the cylinder is analogous to the open 
flow near a transition between two cylindrical ducts (regions 
I and II of Fig. 0. We assume that in region I the flow is 

Uo djcl-c^J^+cd 2 ^) 
given by Eq. O and the cylindrical flow in region II is given " led Y x {^c)Ji(^d) - Y\(jd)Ji{^c) ' 
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In order to appreciate the influence of the inner cylinder disk. Pictures were taken using a 5 megapixel Canon digital 
on the disappearance of cylindrical solutions, figures[2|and[3] camera. To investigate the effect of the inner cylinders on 
show the ratio c/R as a function of the swirl S. The effect of VB, we have used three axial fixed rod of radius d = 1 mm, 
a very slender inner cylinder may be to increase (Fig.|21l or d = 2.5 mm and d = 5 mm. 

n 

decrease (Fig. slightly the critical value S c . However, for For convenient comparison with previous works |j , we 
all the situations considered, the critical value S c increases &st studied the dynamical behavior of usual vortex break- 
as long as d is above a threshold. In this case, the emergence down without axis rod . In Fig. ©the experimental results for 
of VB was transferred to larger values of the swirl param- different aspect ratios can be observed. For H/R = 1 , vortex 
eter, and VB was suppressed in the range between the old breakdown does not appear . For H / R= 1>5> vortex break . 
and new critical values of S (see Fig.0. These results show down take place for Re = 94Q whh further increment of 
that the model predicts the suppressing effect of the slender Reyno id s number, the recirculation bubble becomes oscilla- 
cylinders. In the next section, we present the experimental tory (Re = 1732)> and for Re = ^ h disappears . These re- 
resul , s „ f , h a,„e re pe, t „ m e dt „ studyth ee tf ec totlh e r „ ds sulls agree wilh classical work of Escudier (3). For H/R= 2 
in confined flows. and 2- 5 j fi rst one yB, and subsequently two VB, are gener- 
ated when the Reynolds number is increased. It is possible 
to see in figure the interior detailed structure of the recir- 
culation bubble. 



3 Experimental Setup and Results 



The experimental setup consists of an acrylic cylindrical con- We consider the situation in which an axial fixed rod of 
tainer of inner radius R = 40 mm and a rotating top disk at radius d is introduced. For d = 1 mm, the changes in the 
a variable height H rotating with angular velocity Q. (fig- flow are very small in comparison with the situation without 
ure[5}- The fluid used was water dissolutions of glycerin at rod. However, for d = 2.5 mm and d = 5 mm, we observe 
60% in mass, with V = 1 x 10~ 5 m 2 /s. Temperature was kept that the critical Reynolds number to produce vortex break- 
constant at 20°C. The Reynolds number corresponding to down is increased. Experimental results are summarized in 
the rotating top wall, Re = £2R 2 /v, varies between 600 and Tabled In Figs. [8] and [9] the critical Reynolds number for 
2600, with a 1 % error. Four different aspect ratios H/R were the appearance of the first bubble as a function of the aspect 
used: 1, 1.5, 2 and 2.5. The visualization system consists of ratio and the radius of the rod respectively is shown. We note 
a vertical sheet of light with 2 mm of thickness, generated that for values of Re larger that Re c , the size of the bubble 
by two slide projectors. Small quantities of fluorescein were and its dynamics is still affected by the presence of the inner 
injected into the flow through a small hole in the bottom cylinder. For example, Figure Q0] shows vortex breakdown 
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for H/R = 2.5 and Re = 2260. We can observe that without was d/R = 0.04, which is similar to the small radius ratio 

an axial rod, there are two oscillating bubbles (Fig. UOfa). value we considered, that is d/R = 0.025. Significant results, 

For d = 2.5 mm, there are two VB, but, in this case, they are in our experiments, are obtained above d/R = 0.0625. 

steady (Fig.[l0j)) and for d = 5 mm, only one VB appears. It is also interesting to compare our results with those 

In addition, the size of the bubble corresponding to the first obtained by Mullin et al JlJ). These authors reported that 

VB was clearly decreased with the presence of the rods. they did not obtain noticeable changes in the emergence of 

VB with the addition of a inner cylinder with ratio d/R = 

0.1. However, in our experiments we observed appreciable 

4 Discussion 

changes in the flow, for values of the ratio d/R = 0.0625 and 



0, 



From the tables and figures shown above it can be concluded d l R = 125 - Looking carefully at Figures 2 and 3 of ( 12), we 

that the usage of the inner rods at rest, may be employed observe sli S ht differences which, according to the authors, 

to control the vortex onset. Comparing Figs. 0] and E we are within the experimental errors. For this reason, we took 

observe that the experimental results about the effect of the specially careful measurements in these cases, 
cylinders on the first bubble formation are in agreement with 

the prediction of the simple model, assuming the hypothesis 

5 Summary and Conclusions 

that the swirl parameter increases as Re increases. This de- 
pendence of S on Re is suggested by the fact that both quan- In this work, we present a method for controlling the onset 
tities are proportional to the angular momentum of the flow, of VB. It consists basically of the addition of a small cylin- 
Recently, Husain et al proposed arguments that support this der at rest in the axis of the cylindrical container. The ex- 
hypothesis dlli . periments we performed show that this procedure increases 
The changes of the critical Reynolds number due to the the critical Reynolds number for the emergence of VB, and 
rods presence are clearly noticeable, as follows from table 1 . consequently suppresses the onset of VB in a certain range 
We also observe that in the presence of two bubbles, larger of values of Re. 

values of the Reynolds number are needed for the devel- The control technique proposed here is simpler than other 

opment of the second bubble in comparison with the situ- previously proposed in the literature, since it does not re- 

ation without rod. It is worth noting that our results are in quire additional auxiliary devices. It is worth noting that 

agreement with those of Hussain et al dill , in so far as that our method, unlike the approach proposed by Husain et al 

the very slender rod at rest does not introduce significant dlli . does not imply the addition of swirl near the axis. The 

changes in the flow. The radius ratio used by these authors simplicity of a method is in general an interesting feature, 
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becoming more feasible to be used in engineering devices. 
Moreover, the required modification of the duct is relatively 
small. The volume ratio (cylinder to row) is V1/V2 ~4x 
10~ 3 , while the decrease of the critical Reynolds number is 
about 10% . So that the shift of the critical Reynolds number 
is 20 times larger that the percent modification of the volume 
of the cylinder, showing the effectiveness of the method. 

The effect of the rods is consistent with the results of a 
simple theoretical model of VB based on the failure of the 
quasi-cylindrical approximation, which assumes that S in- 
creases with increasing Re. This agreement supports the idea 
of that to explain VB a local description of the flow is ade- 
quate, and the details of the flows far for the localization of 
VB do not play an important role. This point is very appeal- 
ing from a theoretical point of view, since it implies that to 
study VB in closed flows we must not consider all the rather 
complex velocity fields, but that it is enough to consider only 
the local field near the axis of the cylinder. 

We acknowledge financial support from the Programa de 
Desarrollo de Ciencias Basicas (PEDECIBA, Uruguay) and 
Grants FCE 9028 and PDT54/037 (Conicyt, Uruguay). 
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d=l mm 


d=2.5 mm 
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Fig. 1 (a) Streamlines of the enclosed flow prior to the development of 
VB (taken from Ref. [Escudier 1981 ]). The streamlines are divergent 
near the point where VB develops, suggesting that locally the flow may 
be described as the flow inside a duct that experiment an expansion as 
it is shown in (b). 
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Fig. 2 Theoretical results for the dimensionless core size, c/R, as a 
function of the swirl, S, and different radius of the fixed rod, (a) d = 0, 
(b) d = 0.2, (c) d = 0.4, (d) d = 0.6 and (e) , d = 0.8. Other parameter 
values: R = 1, c = 0.25 and/? = 1.2 



10 



Rcmin? disc 



]{ 



0.U- 





1 J n 



















Fig. 5 Experimental setup. Closed cylinder, radius R and height H, 
with a rotating top wall. Along the axis of the cylinder a fixed rod of 
radius d is located. 
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Fig. 3 As figure|2|but with R = 1, c = 0.4 and R = 1 .2, for (a) d = 0, 
(b) = 0.2, (c) d = 0.4, (d) d = 0.6. 
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Fig. 6 Critical Reynolds number corresponding to the appearance of 
a single VB (circles), double VB (squares), oscillatory bubbles (dia- 
mond) and disappearance of VB (stars) for the case without fixed rod. 




Fig. 4 Critical value of the swirl parameter, 5 C , for the disappearance 
of cylindrical solutions as a function of d/R, for the case Rq = 1, cq = 
0.4 and R= 1.2 




Fig. 7 Flow visualization showing one and two vortex breakdown, 
without axial rod. (a) Re = 1300, H/R = 2. (b) Re = 1756, H/R=2.5. 



Fig. 8 Critical Reynolds numbers corresponding to the appearance of 
the first bubble as a function of the aspect ratio H/R. Without axial 
rod (circles) and with axial rod; d = 2.5 mm (squares) and d = 5 mm 
(triangles). 
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Fig. 9 Critical Reynolds number as a function d/R for the appearance 
of single VB and different aspect ratios, H/R = 2 (squares) and H/R = 
2.5 (circles). 




Fig. 10 Vortex breakdown for H/R = 2.5 and Re=2260. (a) without 
axial rod; (b) for d = 5 mm; (c) for d = 10 mm. 



